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Abstract. We give the first algorithmic study of a class of "covering tour" problems related 
to the geometric Traveling Salesman Problem: Find a polygonal tour for a cutter so that it sweeps 
out a specified region ("pocket"), in order to minimize a cost that depends mainly on the number of 
turns. These problems arise naturally in manufacturing applications of computational geometry to 
automatic tool path generation and automatic inspection systems, as well as arc routing ("postman") 
problems with turn penalties. We prove the NP-completeness of minimum-turn milling and give 
efficient approximation algorithms for several natural versions of the problem, including a polynomial- 
time approximation scheme based on a novel adaptation of the m-guillotine method. 
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1. Introduction. An important algorithmic problem in manufacturing is to 
compute effective paths and tours for covering ( "milling" ) a given region ( "pocket" ) 
with a cutting tool. The objective is to find a path or tour along which to move a 
prescribed cutter in order that the sweep of the cutter covers the region, removing all 
of the material from the pocket, while not "gouging" the material that lies outside 
of the pocket. This covering tour or "lawn mowing" problem 6 and its variants 
arise not only in NC machining applications but also in automatic inspection, spray 
painting/coating operations, robotic exploration, arc routing, and even mathematical 
origami. 

The majority of research on these geometric covering tour problems as well as on 
the underlying arc routing problems in networks has focused on cost functions based on 
the lengths of edges. However, in many actual routing problems, this cost is dominated 
by the cost of switching paths or direction at a junction. A drastic example is given by 
fiber-optical networks, where the time to follow an edge is negligible compared to the 
cost of changing to a different frequency at a router. In the context of NC machining, 
turns represent an important component of the objective function, as the cutter may 
have to be slowed in anticipation of a turn. The number of turns ("link distance") 
also arises naturally as an objective function in robotic exploration (minimum-link 
watchman tours) and in various arc routing problems, such as snow plowing or street 
sweeping with turn penalties. R. Klein |35j has posed the question of minimizing the 
number of turns in polygon exploration problems. 
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In this paper, we address the problem of minimizing the cost of turns in a covering 
tour. This important aspect of the problem has been left unexplored so far in the 
algorithmic community, and the arc routing community has examined only heuristics 
without performance guarantee, or exact algorithms with exponential running time. 
Thus, our study provides important new insights and a better understanding of the 
problems arising from turn cost. We present several new results: 

(1) We prove that the covering tour problem with turn costs is NP-complete, 
even if the objective is purely to minimize the number of turns, the pocket is 
orthogonal (rectilinear), and the cutter must move axis-parallel. The hard- 
ness of the problem is not apparent, as our problem seemingly bears a close 
resemblance to the polynomially solvable Chinese Postman Problem; see the 
discussion below. 

(2) We provide a variety of constant-factor approximation algorithms that ef- 
ficiently compute covering tours that are nearly optimal with respect to 
turn costs in various versions of the problem. While getting some 0{1)- 
approximation is not difficult for most problems in this class, through a careful 
study of the structure of the problem, we have developed tools and techniques 
that enable significantly stronger approximation results. 

One of our main results is a 3.75-approximation for minimum-turn axis- 
parallel tours for a unit square cutter that covers an integral orthogonal poly- 
gon, possibly with holes. Another main result gives a 4/3- approximation for 
minimum-turn tours in a "thin" pocket, as arises in the arc routing version 
of our problem. 

Table n~n summarizes our results. The term "coverage" indicates the number 
of times a point is visited, which is of interest in several practical applications. 
This parameter also provides an upper bound on the total length. 

(3) We devise a polynomial-time approximation scheme (PTAS) for the covering 
tour problem in which the cost is given as a weighted combination of length 
and number of turns, i.e., the Euclidean length plus a constant C times the 
number of turns. For an integral orthogonal polygon with h holes and A'^ 
pixels, the running time is 

20(h) . 7vO(C)^ rp^^ PTAS involves an extension of 
the m-guillotine method, which has previously been applied to obtain PTAS's 
in problems involving only length ^38". 
We should stress that our paper focuses on the graph-theoretic and algorithmic 
aspects of the turn-cost problem; we make no claims of immediate applicability of our 
methods for NC machining. 

Related Work. In the CAD community, there is a vast literature on the subject 
of automatic tool-path generation; we refer the reader to Held [27] for a survey and 
for applications of computational geometry to the problem. The algorithmic study of 
the problem has focused on the problem of minimizing the length of a milling tour: 
Arkin, Fekete, and Mitchell E| show that the problem is NP-hard for the case 
where the mower is a square. Constant-factor approximation algorithms are given in 
[51I^I3L)|. with the current best factor being a 2.5-approximation for min-length milling 
(11/5-approximation for orthogonal simple polygons). For the closely related lawn 
mowing problem (also known as the "traveling cameraman problem" [301), in which 
the covering tour is not constrained to stay within P, the best current approximation 
factor is 3 -I- e (utilizing PTAS results for TSP). Also closely related is the watchman 
route problem with limited visibility (or "d-sweeper problem"); Ntafos 42 provides 
a 4/3-approximation, and Arkin, Fekete, and Mitchell 6 improve this factor to 6/5. 
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Table 1.1 



Approximateort factors achieved by our (polynomial-time) algorithms. Rows marked "APX" 
give approximation factors for the minimum-turn cycle cover ("Cycle cover"), the minimum-turn 
covering tour ("Tour"), and the simultaneous approximation of length of a covering tour ( "Length"). 
The row marked "Max Cover" indicates the maximum number of times a point is visited. The 
parameter 5 denotes the maximum degree in the underlying graph, while the parameter p is the 
maximum number of directions in the graph. The two rows for running time refer to an "explicit" 
description of input pixels and output and a more compact "implicit" encoding of pixels and output. 
(See Section 2 for more detailed definitions.) 



The problem is also closely related to the Hamiltonicity problem in grid graphs; the 
results of |44j suggest that in simple polygons, minimum-length milling may in fact 
have a polynomial-time algorithm. 

Covering tour problems are related to watchman route problems in polygons, 
which have received considerable attention in terms of both exact algorithms (for 
the simple polygon case) and approximation algorithms (in general); see (SHI for a 
relatively recent survey. Most relevant to our problem is the prior work on minimum- 
link watchman tours: see [21 El EI for hardness and approximation results, and |14[I36 | 
for combinatorial bounds. However, in these problems the watchman is assumed to 
see arbitrarily far, making them distinct from our tour cover problems. 

Other algorithmic results on milling include a study of multiple tool milling by 
Arya, Cheng, and Mount 0, which gives an approximation algorithm for minimum- 
length tours that use different size cutters, and the paper of Arkin, Held, and Smith [J], 
which examines the problem of minimizing the number of retractions for "zig-zag" 
machining without "re-milling" , showing that the problem is NP-complete and giving 
an 0(l)-approximation algorithm. 

Geometric tour problems with turn costs have been studied by Aggarwal et al. j^, 
who study the angular-metric TSP. The objective is to compute a tour on a set of 
points, such that the sum of the direction changes at vertices is minimized: For 
any vertex Vi with incoming edge Vi-iVi and outgoing edge WiWi+i, the change of 
direction is given by the absolute value of the angle between Vi-iVi and ViVi^i. The 
problem turns out to be NP-hard, and an O(logn) approximation is given. Fekete |2U | 
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and Fekete and Woeginger |i2T] have studied a variety of angle-restricted tour (ART) 
problems. Covering problems of a different nature have been studied by Demaine et 
al. |15) . who considered algorithmic issues of origami. 

In the operations research literature, there has been an extensive study of arc 
routing problems, which arise in snow removal, street cleaning, road gritting, trash 
collection, meter reading, mail delivery, etc.; see the surveys of 1181 Arc 
routing with turn costs has had considerable attention, as it enables a more accurate 
modeling of the true routing costs in many situations. Most recently, Clossey, Laporte, 
and Soriano |13| present six heuristic methods of attacking arc routing with turn 
penalties, without resorting to the usual transformation to a TSP problem; however, 
their results are purely based on experiments and provide no provable performance 
guarantees. The directed postman problem in graphs with turn penalties has been 
studied by Benavent and Soler , who prove the problem to be (strongly) NP-hard 
and provide heuristics (without performance guarantees) and computational results. 
(See also Fernandez's thesis |22| and for computational experience with worst-case 
exponential-time exact methods.) 

Our covering tour problem is related to the Chinese Postman Problem (CPP), 
which can be solved exactly in polynomial time for "purely" undirected or purely 
directed graphs. However, the turn-weighted CPP is readily seen to be NP-complete: 
Hamiltonian cycle in line graphs is NP-complete (contrary to what is reported in |25| : 
see page 246 of West ^^), implying that TSP in Hue graphs is also NP-complete. 
The CPP on graph G with turn costs at nodes (and zero costs on edges) is equivalent 
to TSP on the corresponding line graph, C[G), where the cost of an edge in £(G) 
is given by the corresponding turn cost in G. Thus, the turn- weighted CPP is also 
NP-complete. 

2. Preliminaries. This section formally defines the problems at hand, and var- 
ious special cases of interest. 

Problem Definitions. The general geometric milling problem is to find a closed 
curve (not necessarily simple) whose Minkowski sum with a given tool (cutter) is 
precisely a given region P bounded by n edges. In the context of numerically controlled 
(NC) machines, this region is usually called a pocket. Subject to this constraint, we 
may wish to optimize a variety of objective functions, such as the length of the tour, 
or the number of turns in the tour. We call these problems minimum-length and 
minimum-turn milling, respectively. While the latter problem is the main focus of 
this paper, we are also interested in bicriteria versions of the problem in which both 
length and number of turns must be small; we also consider the scenario in which the 
objective function is given by a linear combination of turn cost and distance traveled 
(see Section 

In addition to choices in the objective function, the problem version depends on 
the constraints on the tour. The most general case arises when considering a tour that 
has to visit a discrete set of vertices, connected by a set of edges, with a specified turn 
cost at each vertex to change from one edge to the next. More precisely, at each vertex, 
the tour has the choice of (0) going "straight" , if there is one "collinear" edge with 
the one currently used (costing no turn), (1) turning onto another, non-coUinear edge 
(costing one turn), or (2) "U-turning" back onto the source edge (costing two turns). 
Which pairs of the d edges incident to a degree-d vertex are considered "collinear" 
is specified by a matching in the complete graph Kd'. three incident edges cannot 
be "collinear" . We call this graph-theoretic abstraction the discrete milling problem^ 
indicating the close relationship to other graph-theoretic tour optimization problem. 
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We are able to give a number of approximation algorithms for discrete milling that 
depend on some graph parameters: S denotes the maximum degree of a vertex and 
p denotes the maximum number of distinct "directions" coming together at a vertex. 
(For example, for graphs arising from c?-dimensional grids, these values are bounded 
by 2c? and d, respectively.) 

A special case of discrete milling arises when dealing with "thin" structures in two- 
or three-dimensional space, where the task is to travel all of a given set of "channels" , 
which are connected at vertices. This resembles a CPP, in that it requires us to travel 
a given set of edges; however, in addition to the edge cost, there is a cost at the 
vertices when moving from one edge to the next. For this scenario, we are able to 
describe approximation factors that are independent of other graph parameters. 

More geometric problems arise when considering the milling of a polygonal region 
P. In the orthogonal milling problem, the region P is an orthogonal polygonal domain 
(with holes) and the tool is an (axis-parallel) unit-square cutter constrained to axis- 
parallel motion, with edges of the tour alternating between horizontal and vertical. 
All turns are orthogonal; 90° turns incur a cost of 1, while a "U-turn" has a cost of 2. 
In the integral orthogonal case, all coordinates of boundary edges are integers, so the 
region can be considered to be the (connected) union of TV pixels, i.e., axis-parallel 
unit squares with integer vertices. Note that in general, N may not be bounded by 
a polynomial in n. Instead of dealing directly with a geometric milling problem, we 
often find it helpful to consider a more combinatorial problem, and then adapt the 
solution back to the geometric problem. In particular, for integral orthogonal milling, 
we may assume that an optimal tour can be assumed to have its vertex coordinates of 
the form k + ^ for integral k. Then, milling in an integral orthogonal polygon (with 
holes) is equivalent to finding a tour of all the vertices ( "pixels" ) of a grid graph; see 
Figure O 



I 




Figure 2.1. An instance of the integral orthogonal milling problem (left) and the grid graph 
model (right). 

An interesting special case of integral orthogonal milling is the thin orthogonal 
milling problem, in which the region does not contain a 2x2 square of pixels. This 
is also closely related to discrete milling, as we can think of edges embedded into 
the planar grid, such that vertices and channels are well separated. This problem of 
finding a tour with minimum turn cost for this class of graphs is still NP-complete, 
even for a subclass for which the corresponding problem of minimizing total distance is 
trivial; this highlights the particular difficulty of dealing with turn cost. On the other 
hand, thin orthogonal milling allows for particularly fast and efficient approximation 
algorithms. 
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Other Issues. It should be stressed that using turn cost instead of (or in addition 
to) edge length changes several characteristics of distances. One fundamental problem 
is illustrated by the example in Figure the triangle inequality does not have to 
hold when using turn cost. This implies that many classical algorithmic approaches 
for graphs with nonnegative edge weights (such as using optimal 2-factors or the 
Christofides method for the TSP) cannot be applied without developing additional 
tools. 



a 



• 1 

b 

Figure 2.2. The triangle inequality may not hold when using turn cost as distance measure: 
d(a, c) = 3 > 2 = d{a, b) + d(h, c) 



In fact, in the presence of turn costs we distinguish between the terms 2-factor, i.e., 
a set of edges, such that every vertex is incident to two of them, and cycle cover, i.e., 
a set of cycles, such that every vertex is covered. While the terms are interchangeable 
when referring to the set of edges that they constitute, we make a distinction between 
their respective costs: a "2-factor" has a cost consisting of the sum of edge costs, but 
does not necessarily account for the turn cost between it two incident edges, while 
the cost of a "cycle cover" includes also the turn costs at vertices. 

It is often useful in designing approximation algorithms for optimal tours to begin 
with the problem of computing an optimal cycle cover, minimizing the total number 
of turns in a set of cycles that covers P. Specifically, we can decompose the problem 
of finding an optimal (minimum-turn) tour into two tasks: finding an optimal cycle 
cover, and merging the components. Of course, these two processes may influence 
each other: there may be several optimal cycle covers, some of which are easier to 
merge than others. (In particular, we say that a cycle cover is connected, if the graph 
induced by the set of cycles and their intersections is connected.) As we will show, 
even the problem of optimally merging a connected cycle cover is NP-complete. This 
is in contrast to minimum-length milling, where an optimal connected cycle cover can 
trivially be converted into an optimal tour that has the same cost. 

Algorithms whose running time is polynomial in the explicit encoding size (pixel 
count) are ■pseudo-polynomial. Algorithms whose running time is polynomial in the 
implicit encoding size are polynomial. This distinction becomes an important issue 
when considering different ways to encode input and output; e.g., a large set of pix- 
els forming an a x 6 rectangle can be described in space 0(loga -I- log 6) by simply 
describing the bounding edges, instead of listing all ab individual pixels. In integral 
orthogonal milling, one might think that it is most natural to encode the grid graph 
with vertices, because the tour will be embedded on this graph and will, in general, 
have complexity proportional to the number N of pixels. But the input to any ge- 
ometric milling problem has a natural encoding by specifying only the n vertices of 
the polygon P. In particular, long edges are encoded in binary (or with one real 
number, depending on the model) instead of unary. It is possible to get a running 
time depending only on this size, but of course we need to allow for the output to 
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be encoded implicitly. That is, we cannot explicitly encode each vertex of the tour, 
because there are too many (the number can be arbitrarily large even for a succinctly 
encodable rectangle). Instead, we encode an abstract description of the tour that is 
easily decoded. 

Finally, we mention that many of our results carry over from the tour (or cycle) 
version to the path version, in which the cutter need not return to its original position. 
In this paper, we omit the straightforward changes necessary to compute optimal 
paths. A similar adjustment can be made for the related case of lawn mowing, in 
which the sweep of the cutter is allowed to go outside P during its motion. Clearly, 
our techniques are also useful for scenarios of this type. 

3. NP-Completeness. Arkin, Fekete, and Mitchell [HI have proved that the 
problem of optimizing the length of a milling tour is NP-hard. Their proof is based 
on the well-known hardness of deciding whether a grid graph has a Hamiltonian cycle 
|29[ IST| . This result implies that it is NP-hard to find a tour of minimum total length 
that visits all vertices. If, on the other hand, we are given a connected cycle cover of 
a graph that has minimum total length, then it is trivial to convert it into a tour of 
the same length by merging the cycles into one tour. 

In this section we show that if the quality of a tour is measured by counting turns, 
then even this last step of turning an optimal connected cycle cover into an optimal 
tour is NP-complete. Thus we prove that it is NP-hard to find a milling tour that 
optimizes the number of turns for a polygon with holes. 

Theorem 3.1. Minimum-turn milling is NP-complete, even when we are re- 
stricted to the orthogonal thin case, and are already provided with an optimal connected 
cycle cover. 

Because thin orthogonal milling is a special case of thin milling as well as orthog- 
onal milling, and it is easy to convert an instance of thin orthogonal milling into an 
instance of integral orthogonal milling, we have 

Corollary 3.2. Discrete milling, orthogonal milling, and integral orthogonal 
milling are NP-complete. 

Proof. [Theorem 13. 1| Our reduction proceeds in two steps. First we show that 
the problem Hamiltonicity of Unit Segment Intersection Graphs (Husig) 
of deciding the Hamiltonicity of intersection graphs of axis-parallel unit segments is 
hard. To see this, we use the NP-hardness of deciding Hamiltonicity of grid graphs 
f |29[l31p and argue that any grid graph can be represented in this form; see Figure l^TI 

Consider a set of integer grid points that induce a grid graph G. Note that G 
is bipartite, because one can 2-color the nodes by coloring a grid point {x, y) black 
(rcsp., white) if a; -I- y is odd (resp., even). After rotating the point set by 7r/4, the 
coordinate of each point is an integer multiple of 1/ V2. Scaling down the resulting 
arrangement by a factor of 3/\/2 results in an arrangement in which the coordinate 
of each point is an integer multiple of 1/3, and the shortest distance between two 
points of the same color class is 2/3. For the resulting set of points pi — {xi,yi), let 
p'i = Pi + be given as the set obtained by "perturbations" £i that are small 

and all distinct. Then, represent each "white" vertex by a horizontal unit segment 
centered at p[, and each "black" vertex by a vertical unit segment centered at p'^. 
Now it is easy to see that the resulting unit segment intersection graph is precisely 
the original grid graph G. 

In a second step, we show that the problem HusiG reduces to the problem of 
milling with turn costs. The outline of our argument is illustrated in Figure IX^ 

Consider a unit segment intersection graph G, given by a set of axis-parallel unit 
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(a) (b) 

Figure 3.1. (a) A grid graph G. (b) A representation of G as an intersection graph of axis- 
parallel unit segments. 



segments, as shown in Figure rOT a). Figure |^2Ib) shows the corresponding graph, 
with a Hamiltonian cycle indicated in bold. Without loss of generality, we may assume 
that G is connected. Let s be the number of nodes of G. 

As shown in Figure rOT c) . we replace each line segment by a cycle of four thin 
axis-parallel corridors. This results in a connected polygonal region P having 4s 
convex corners. Clearly, any cycle cover or tour cover of P must have at least 4s 
turns; by using a cycle for each set of four corridors representing a strip, we get a 
cycle cover C with 4s turns. Therefore, C is an optimal cycle cover, and it is connected, 
because G is connected. 

Now assume that G has a Hamiltonian cycle. It is easy to see (Figure I3.2r f)') 
that this cycle can be used to construct a milling tour of P with a total of 5s turns: 
Each time the Hamiltonian cycle moves from one vertex vi of the grid graph to the 
next vertex Vj , the milling tour moves from the cycle Ci representing Vi to the cycle 
Cj representing Wj, at an additional cost of 1 turn for each of the s edges in the 
Hamiltonian cycle. 

Assume conversely that there is a milling tour T with at most 5s turns. We refer 
to turns at the corners of 4-cycles as convex turns. The other turns are called crossing 
turns. 

As noted above, the convex corners of P require at least 4s convex turns. Con- 
sider the sequence of turns ii, . . . , in T . By construction, the longest contiguous 
subsequence of convex turns contains at most four different convex turns. (More pre- 
cisely, we can only have such a subsequence with four different convex corners, if these 
four corners belong to the same 4-cycle representing a unit segment.) Furthermore, 
we need at least one additional crossing turn at an interior crossing of two corridors to 
get from one convex corner to another convex corner not on the same 4-cycle. (More 
precisely, one crossing turn is sufficient only if the two connected convex corners be- 
long to 4-cycles representing intersecting unit segments.) Therefore, we need at least 
c crossing turns if we have at least c contiguous subsequences as described above. 
This means that c > s; hence, c = s by the assumption on the number of turns on T. 
Because the c crossing turns correspond to a closed roundtrip in G that visits all s 
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successor 



(e) (f) 

Figure 3.2. Thin orthogonal milling with turn cost is NP-hard: (a) a set of s = 8 axis- 
parallel unit segments, denoted by si, . . . , sg; (b) the corresponding intersection graph G, with the 
Hamiltonian cycle vi,V2,V3,Vfi,VY,Vfi,v$,V4 shown in bold; (c) representing G by a connected region 
consisting of is corridors; (d) a drawing of the graph induced by the instance of thin orthogonal 
milling, with the s = 8 rectangular cycles Ci, . . . , Cg; (e) a milling tour with bs turns corresponding 
to the Hamiltonian cycle in G; (f) milling the four corridors of a cycle using five turns. 

vertices, this implies that we have a Hamihonian cycle, concluding the proof. □ 

4. Approximation Tools. There are three main tools that we use to develop 
approximation algorithms: computing optimal cycle covers for milling the "boundary" 
of P ('Section |4.1|l . converting cycle covers into tours f Section |4. 2(1 . and using optimal 
(or nearly-optimal) "strip covers" fSection l4.3|l . In this section, our description mostly 
focuses on orthogonal milling; however, we will see in the following Section 5.1 how 
some of our tools can also be applied to the general case of discrete milling. 
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4.1. Boundary Cycle Covers. We consider first the problem of finding a 
minimum-turn cycle cover for covering a certain subset, P, of P that is along its 
boundary. This will turn out to be a useful tool for approximation algorithms. Specif- 
ically, in orthogonal milling we define the set P of boundary pixels to consist of pixels 
that have at least one of their four edges on a boundary edge of the polygon; i.e., in 
the grid graph that describes adjacency of pixels, these are pixels of degree at most 
3. Let N-p be the number of boundary pixels. A boundary cycle cover is a collection 
of cycles that visit all boundary pixels. 

We define an auxiliary structure, Gp = {V-p,E-p), which is a complete weighted 
graph on 2N^ vertices; for ease of description, we will refer to Gp as a set of points 
and paths between them. This will allow us to map boundary cycle covers in P to 
matchings of corresponding turn cost in G-p. For this purpose, map each pixel pi € P 



to two vertices in V-p, v. 



(0) 



and 



.,(1) 



For each boundary pixel pi , this pair represents an 



orientation that is attained by a cutter when visiting pi. Depending on the boundary 
structure oi pi, there are four different cases; refer to Figure ^T] 
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(iv) 



Figure 4.1. Representing a boundary pixel pi by a pair of vertices i;^"' and v^^'^ . 



(i) One edge of pi is a boundary edge of the polygon. 

(ii) Two opposite edges of pi are boundary edges of the polygon. 

(iii) Two adjacent edges of pi are boundary edges of the polygon. 

(iv) Three edges of Pi are boundary edges of the polygon. 



For easier description, we refer to the vertices vf''^ and v^''' as points embedded 
within Pi, as shown in Figure im Furthermore, we add a mandatory path rui between 
vf^ and v''^'\ represented by a polygonal path with c{mi) = (cases (i) and (ii)), 
c{mi) = 1 (case (iii)) or c{mi) = 2 turns (case (iv)), as shown in the figure. This 
path maps the contour of P at pi, and it represents orientations that a cutter has to 
attain when visiting pixel p^. Note that traveling from v^'^'^ to Vj^ along induces 
a heading (5,[''' •* when leaving v\'^^ and a heading 6^^ when arriving at 
Note that S^'^' is opposite to Si'^'^K 

Now we add a set of optional paths, representing the weighted edges Ep of the 
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complete graph G-p. For an example, refer to Figure 1^?^ For any pair of vertices v, 



(h) 



and Vj'^'^ , let d{v'^'^' , Vj'^' ) be the minimum number of turns necessary when traveling 
from Wj-'*'' with heading S^'^''^^ to vf' with heading (jj''' Note that d{v^'^\ v^'''' ) = 



.W Jk) 



d{v 



(fc) [h) 



), as any shortest path can be traveled in the opposite direction. Using 



a Dijkstra-like approach, we can compute these distances from one boundary pixel 
to all other boundary pixels in time O {N-p log N-p); see the overview in PH] or the 
paper The overall time of 0{N^logN-p) for computing all these link distances 
is dominated by the following step: In time 0{N^) find a minimum-weight 

perfect matching in the complete weighted graph Gp. 



H---J----J- 



5(0,-) 



,(0) 



Figure 4.2. The cost of traveling between two pixels: d{v^^\v^_^^) = d{v^^\v^^^^) = 2, while 



d{v\ 



(0) 



Now it is not hard to see the following. 

Lemma 4.1. Any boundary cycle cover in P with t turns can be mapped to a 
perfect matching in G-p of cost t — '^p.^-pc{mi), and vice versa. 

Proof. Whenever a boundary cycle cover visits a boundary pixel p^, it has to 
perform the turns corresponding to the mandatory path rui. Moreover, moving from 
one pixel pi to the next pixel pj can be mapped to the optional path corresponding 



to the edges [v. 



(h) Jk) 



); clearly, the overall cost is as stated. 



Conversely, it is straightforward to see that the combination of a perfect matching 
in Gp and the mandatory paths yields a boundary cycle cover of corresponding cost. 



Using the algorithms described above, we also obtain the following: 

Theorem 4.2. Given the set of N-p boundary pixels, a minimum-turn boundary 

cycle cover can be computed in time 0{N-p), the time it takes to compute a perfect 

matching in Gp. 

If the set of pixels is not given in unary, but implicitly as the pixels contained in 
a region with n edges, the above complexity is insufficient. However, we can use local 
modifications to argue the following tool for speeding up the search for an optimal 
perfect matching. 

Lemma 4.3. Letpi and pj be neighboring boundary pixels that are adjacent to the 
same boundary edge, so d{vf^\v'^^^) — for an appropriate choice of h and k. Then 
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there is an optimal matching containing (v^ \vj ) . 

Proof. This foUows by a simple exchange argument. See Figure ^3] Suppose 
two adjacent pixels pi and pj along the same boundary edge are not matched to each 
other, let w^''-* be the vertex such that Si{h, — ) is heading for pj, and let Wj*^^ be the 
vertex such that 6j{k, —) is heading for pi. Furthermore suppose that v^!^ ■* is matched 

to u,-'''' and u^f is matched to v'^''\ Then we can match v^!^ to '' and to v^''^ 
without changing the cost of the matching. □ 
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Figure 4.3. By performing local modifications, an optimal cycle cover can be assumed to cover 
each collinear piece of the boundary in one connected strip. 



This allows us to obtain a strongly polynomial version of the matching algorithm 
of Theorem lO 

Theorem 4.4. A minimum-turn boundary cycle cover can he computed in time 

Proof. By applying Lemma 14.31 repeatedly, we get 0{n) connected boundary 
strips, consisting of sets of collinear boundary pixels. These can be determined effi- 
ciently by computing offsets of the boundary edges. This leaves only 0{n) endpoints 
of such strips to be matched, resulting in the claimed complexity. □ 

Note that the validity of this argument is not restricted to the integral orthogonal 
case, but remains valid even for orthogonal regions with arbitrary boundary edges. 

Remark. The definition of the "boundary" pixels P used here does not include 
all pixels that touch the boundary of P in a diagonal fashion; in particular, it omits 
the "reflex pixels" that share a corner, but no edge, with the boundary of P. It seems 
difficult to require that the cycle cover mill reflex pixels, because Lemma fl.Hl does not 
extend to this case, and an optimal cycle cover of the boundary (as defined above) 
may have fewer turns than an optimal cycle cover that mills the boundary P plus the 
reflex pixels; see Figure ^31 

4.2. Merging Cycles. It is often easier to find a minimum-turn cycle cover (or 
constant-factor approximation thereof) than to find a minimum-turn tour. We show 
that an exact or approximate minimum-turn cycle cover implies an approximation for 
a minimum-turn tour. 

We concentrate on the integral orthogonal case. First we define a few terms 
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Figure 4.4. Optimally covering all pixels that have an edge against the boundary can leave 
reflex pixels uncovered. 

precisely. Two pixels are adjacent if the distance between their centers is 1. Two 
cycles Ti, T2 are intersecting iff Ti n T2 7^ 0. Two cycles are called touching if and 
only if they are not intersecting and there exist pixels pi € Ti, p2 S T2 such that pi 
and p2 are adjacent. 

Lemma 4.5. Let Pi and P2 be two cycles, with ti andt2 turns, respectively, and 
let p be a pixel that is contained in both cycles. Then there is a cycle milling the union 
of pixels milled by Pi and P2 and having at most ti + 12 + 2 turns. This cycle can be 
found in time linear in the number of its turns. 

Proof. Let the neighbors of p in Pi be ai, 02 and those of p in P2 be 61, 62. Connect 
ai via p to bi and 02 via p to 62 to get the required tour. The two connections may 
add at most a turn each. Hence the resulting tour can be of size at most <i + <2 + 2. □ 

Lemma 4.6. Given two touching cycles Ti, T2 with ti, t2 turns respectively, then 
there is a tour T with at most ti + t2 + 2 turns that mills the union of pixels milled 
by Ti, T2. 

Proof. Because Ti, T2 are touching, Ti n 72 = and there exist adjacent pixels 
pi e Ti and p2 € 12. Without loss of generality assume that p2 is a leftmost such 
pixel, and below pi. Due to these constraints, T2 can enter/exit p2 from only two 
sides. Hence there are only three ways in which T2 can visit p2. These are shown in 
Figure H31 For all three ways we show in Figure how to cut and extend tour T2 
without adding any extra turns, to get a path P2 starting and ending at pixel pi. Cut 
Ti at pi to get a path Pi. By possibly adding two turns, we can merge the two paths 
into one tour. □ 

With the help of these lemmas, we deduce the following: 

Theorem 4.7. A cycle cover with t turns can be converted into a tour with at 
most t + 2(c — 1) turns, where c is the number of cycles. 

Proof. We prove this theorem by induction on the number of tours, c, in the cycle 
cover. The theorem is trivially true for c = 1. For any other c, choose any c — 1 
cycles with a total of t' turns and find a tour T' that covers those c — 1 cycles; by 
induction, it has t' + 2(c — 1) turns. Let the remaining cycle, R, have r turns. Thus 
t = t' + r. Because the polygon is connected, the set of pixels milled by R and T' must 
be connected. Hence either T' and R are intersecting or touching. By Lemmas 14.51 
and 14. til we can merge R and T' into a single tour T with at most t' + 2(c— l) + r + 2 
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T P 



(a) 



(b) 



(c) 



^1 f f 



Figure 4.5. Merging two touching tours: There are three possible ways of tour T2 visiting pixel 
P2 (above). In each case, we can modify T2 into a path that visits pixel pi (below); at a cost of 
possibly one extra turn at each end of the path, we can merge it with tour T\ . 



turns, i.e. t + 2c turns. □ 
Corollary 4.8. A cycle cover of a connected rectilinear polygon with t turns 

can he converted into a single milling tour with at most |t turns. 

Proof. This follows immediately from Theorem 14.71 and the fact that each cycle 

has at least four turns. □ 
Unfortunately, general merging is difficult (as illustrated by the NP-hardness proof 

of Theorem l3.1|) . so we cannot hope to improve these general merging results by more 

than a constant factor. 

4.3. Strip and Star Covers. A key tool for approximation algorithms is a 
covering of the region by a collection of "strips." A strip is a maximal straight 
segment whose Minkowski sum with the tool is contained in the region. A strip cover 
is a collection of strips whose Minkowski sums with the tool cover the entire region. 
A minimum strip cover is a strip cover with the fewest strips. 

Lemma 4.9. The size of a minimum strip cover is a lower bound on the number 
of turns in a cycle cover (or tour) of the region. 

Proof. Any cycle cover induces a strip cover by extending each edge to have 
maximal length. The number of strips in this cover equals the number of turns in the 
cycle cover. □ 

In the discrete milling problem, a related notion is a "rook placement." A rook 
is a marker placed on a pixel, which can attack every pixel to which it is connected 
via a straight axis-parallel path inside the region. A rook placement is a collection of 
rooks no two of which can attack each other. See Figure for an illustration; this 
tool will be used in Theorem 15. 61 based on the following lemma. 

Lemma 4.10. The size of a maximum rook placement is a lower hound on the 
number of turns in a cycle cover ( or tour) for discrete milling. 

Proof. Consider a rook placement and a cycle cover of the region, which must 
in particular cover every rook. Suppose that one of the cycles visits rooks qi, . . . ,qk 
in that order. No two rooks can be connected by a single straight axis-parallel line 
segment, so the cycle must turn inbetween each rook, for a total of at least k turns. 
Because each rook is traversed by at least one cycle, the number of turns (and hence 
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the number of segments in a tour) is at least the number of rooks. □ 
In the integral orthogonal milling problem, the notions of strip cover and rook 

placement are dual and efficient to compute: 

Lemma 4.11. For integral orthogonal milling, a minimum strip cover and a 

maximum rook placement have equal size. For a polygonal region with n edges and N 

pixels they can he computed in time 0{N^'^'^^) or 0{n^'^ logN). 
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Figure 4.6. (Left) An orthogonal region, its subdivision into axis-parallel strips, and a resulting 
greedy rook cover (indicated by black pixels). (Right) An optimal strip cover, and an optimal rook 
cover (indicated by black pixels). 

Proof. For the case of iV e 0{n), the claim follows from Proposition 2.2 in j26) : 
We rephrase the rook-placement problem as a matching problem in a bipartite graph 
G = {Vi,V2, E). Let the vertices in Vi correspond to vertical strips, and the vertices 
in V2 correspond to horizontal strips. An edge e — (u,v) (z E exists if the vertical strip 
corresponding to u and the horizontal strip corresponding to v have a pixel in common 
(i.e., the strips cross). It is easy to see that a maximum-cardinality matching in this 
bipartite graph corresponds to a rook placement: each edge {u,v) in the matching 
corresponds to the unique pixel that vertical strip u and horizontal strip v have in 
common. 

Similarly, observe that the minimum strip-cover problem is equivalent to a min- 
imum vertex-cover problem in the bipartite graph defined above. Each strip in the 
strip cover defines a vertex in the vertex cover. The requirement that each pixel must 
be covered by at least one strip is equivalent to the requirement that each edge of the 
graph must be covered by at least one vertex. 

By the famous Konig-Egervary theorem, the maximum cardinality matching in 
a bipartite graph is equal in size to the minimum vertex cover, and therefore both 
can be solved in time polynomial in the size of the graph; more precisely, this can be 
achieved in time 0{N'^), the time needed for multiplying two N x N matrices, for 
example uj — 2.376; see the paper |2H|j or the survey in chapter 16 of 03], which also 
lists other, more elementary methods. 

To get the claimed running time even for "large" N , using implicit encoding, we 
decompose the region into "thick strips" by conceptually coalescing adjacent horizon- 
tal strips with the same horizontal extent, and similarly for vertical strips. In other 
words, thick strips are bounded by two vertices of the region, and hence there are only 
0{n) of them. We define the same bipartite graph but add a weight to each vertex 
corresponding to the width of the strip (i.e., the number of strips coalesced). Instead 
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of a matching, in which each edge of the graph is either included in the matching or 
not, we now have a multipUcity for each edge, which is the minimum of the weights 
of its two endpoints. The interpretation is that an edge corresponds to a rectangle in 
the region (the intersection of two thick strips) , and the number of rooks that can be 
placed in such a rectangle is at most the minimum of its width and height. 

The weighted-matching problem we consider is that each edge can be included 
in the matching with a multiplicity up to its weight. Furthermore, the sum of the 
included multiplicities of edges incident to a vertex cannot exceed the weight of the 
vertex. A weighted version of the Konig-Egervary theorem states that the minimum- 
weight vertex cover is equal to the maximum- weight matching. (This weighted version 
can be easily proved using the max-flow min-cut theorem.) Both problems can be 
solved in polynomial time using a max-flow algorithm, on a modified graph in which 
a source vertex s is added with edges to all vertices in Vi, of capacity equal to the 
weight of the vertex, and a sink vertex t is added with edges to it from all vertices 
in V2 with capacity equal to the vertex capacity. Edges between Vi and V2 are 
directed from Vi and have capacity equal to the weight of the edge. Currently, the 
best known running time is 0{y/nmlognW) for a bipartite graph with n vertices, m 
edges, and maximum weight W [2'^\ l4Hj . For our purposes, this yields a complexity 
of 0(n2-MogiV). □ 

Note that using weights on the edges is crucial for the correctness of our objective; 
moreover, this has a marked effect on the complexity of the problem: Finding a 
minimum number of axis-parallel rectangles (regardless of their size) that covers an 
integral orthogonal polygon is known to be an NP-complete problem, even for the 
case of polygon without holes '151. 

For general discrete milling, it is possible to approximate an optimal strip cover 
as follows. Greedily place rooks until no more can be placed (i.e., until there is no 
unattackable vertex). This means that every vertex is attackable by some rook, so 
by replacing each rook with all possible strips through that vertex, we obtain a strip 
cover of size p times the number of rooks, where p is the maximum degree of the 
underlying graph. (We call this type of strip cover a star cover.) But each strip in a 
minimum strip cover can only cover a single rook, so this is a p-approximation to the 
minimum strip cover. We have thus proved 

Lemma 4.12. In discrete milling, the number of stars in a greedy star cover is a 
lower bound on the number of strips, and hence serves as a p- approximation algorithm 
for minimum strip covers. Computing a greedy star cover can be done in time 0{N). 

Proof. Loop over the vertices of the underlying graph. Whenever an unmarked 
vertex is found, add it to the list of rooks, and mark it and all vertices attackable by 
it. Now convert each rook into a star as in the proof of Lemma f4.1UI Each edge is 
traversed only once during this process. □ 

5. Approximation Algorithms. We employ four main approaches to building 
approximation algorithms, repeatedly in several settings: 

(i) Star cover + doubling + merging 

The simplest but most generally applicable idea is to cover the region by a 
collection of stars. "Doubling" these stars results in a collection of cycles, 
which can then be merged into a tour using general techniques. 

(ii) Strip cover + doubling + merging 

Tighter bounds can be achieved by covering directly with strips instead of 
stars. Similar doubling and merging steps follow. 

(iii) Strip cover + perfect matching of endpoints + merging 
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The covering of the region is done by the strip cover. To connect these strips 
into cycles, we find a minimum-weight perfect matching on their endpoints. 
This resuhs in a cycle cover, which can be merged into a tour using general 
techniques. 

(iv) Boundary tour + strip cover + perfect matching of odd-degree 
vertices 

Again coverage is by a strip cover, but the connection is done differently. We 
add a tour of the boundary (by merging an optimal boundary cycle cover), 
and attach each strip to this tour on each end. The resulting graph has several 
degree-three vertices, which we fix by adding a minimum-weight matching on 
these vertices. 

5.1. Discrete Milling. As described in the Preliminaries, we consider two sce- 
narios: While general discrete milling focuses on vertices (and thus resembles the 
TSP), thin discrete milling requires traveling a set of edges, making it similar to the 



5.1.1. General Discrete Milling. Our most general approximation algorithm 
for the discrete milling problem runs in linear time. First we take a star cover ac- 
cording to Lemma l4.12l which approximates an optimal strip cover to within a factor 
of p. Then we tour the stars using an efficient method described below. Finally we 
merge these tours using Theorem 14 .71 

We tour each star emanating from a vertex v using the following method — see 
Figure 15.11 Consider a strip s in the star, and suppose its ends are the vertices Wj 
and Uj . A strip having both of its endpoints distinct from v is called a full strip; a 
strip one of whose endpoints is equal to v is called a half strip. Half strips are covered 
by three edges, {v,Uj), {uj,Uj), and {uj,v), making a U-turn at endpoint Uj. (This 
covering is shown for the half strip (y, ui) in Figure lHTT bV) Full strips are covered by 
five edges, {v,Ui), (ui,Ui), {ui,Uj), {uj,Uj), and {uj,v), with U-turns at both ends, 
Ui and Uj. (This covering is shown for the full strip {u^,U2) in Figure [^TT bV) Now 
we have several paths of edges starting and ending at v. By joining their ends we can 
easily merge these paths into a cycle. 



Figure 5.1. (a) A star of degree five around vertex v. {v,ui), [v^u-j,), (i), U4) are half strips, 
{u2,uc,) is a full strip, (b) A covering with three edges for a half strip, and a covering with five edges 
for the full strip. 



CPP. 




The number of turns in this cycle is 3 times the number of half strips, plus 5 times 
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the number of full strips. This is equivalent to the number of distinct directions at v 
plus 2 times the degree of v. (The number of directions at v is equal to the number 
of full strips plus the number of half strips, by definition. The degree of v is equal to 
2 times the number of full strips plus the number of half strips.) Lemma |4 . 1 21 implies 
that the number of stars is a lower bound on the number of turns in a cycle cover of 
the region, proving the following: 

Theorem 5.1. There is an 0{N)-time {25 + p)- approximation for finding a 
minimum-turn cycle cover in discrete milling. Furthermore, the maximum coverage 
of a vertex (i.e., the maximum number of times a vertex is swept) is d, and the cycle 
cover is a 6 -approximation on length. 

Proof. As the star cover, by definition, contains all vertices, the cycle cover 
obtained by traversing the stars also does. As stated above, the number of turns in 
each cycle covering a star is the number of directions at v, plus 2 times the degree of 
V. Summing over all stars, we get the claimed approximation bound. The running 
time follows directly from Lemma 14.121 Deriving the values for maximum coverage 
and overall length is straightforward. □ 

Corollary 5.2. There is a linear-time {2S + p + 2) -approximation for minimum- 
turn discrete milling. Furthermore, the maximum coverage of a vertex is S, and the 
tour is a 6 -approximation on length. 

Proof. We apply Theorem 14. 71 The number of cycles to be merged is the number 
of stars, which by Lemma 14.121 is a lower bound on the number of turns in a tour of 
the region. We pay at most two turns per cycle for the merge. There is no additional 
cost of length due to the merge, as the stars form a connected graph. □ 

5.1.2. Thin Discrete Milling. As described in Scction[51 a more special struc- 
ture arises if the structure to be milled consists of a connected set of "channels" that 
have to be milled. In this case, achieving a strip cover is trivial: 

Lemma 5.3. In thin discrete milling a strip cover can be obtained in linear time 
by merging edges that are collinear at some vertex. 

Using method (ii) described at the beginning of the section, we get the following 
approximation results: 

Theorem 5.4. There is a A- approximation of complexity O(nlogn) for comput- 
ing a minimum-turn cycle cover for a graph with n edges, and a 6 -approximation of 
the same complexity for computing minimum-turn tours. 

Proof. Clearly, the number of strips is a lower bound on the cost of any cycle 
cover or tour. Turning each strip into a cycle with two u-turns, i.e., 4 turns, yields 
a cycle cover within a factor 4 of the optimum. Merging these cycles at a cost of 2 
turns per merge yields a tour within a factor of 6 times the optimum, as each cycle 
has 4 turns. □ 

Clearly, all edges get covered twice (yielding a bound of 2 on the simulataneous 
length approximation) and no vertex gets covered more than 2p times. 

Using the more time-consuming method (iii) , we get better approximation factors 
for the turn cost: 

Theorem 5.5. There is a 1.5 -approximation of complexity 0{n^) for computing 
a minimum-turn cycle cover for a graph with n edges, and a 5. 5 -approximation of the 
same complexity for computing minimum-turn tours. 

Proof. As before, we can compute an optimal strip cover S in linear time. Anal- 
ogous to the approach in Section 14.11 define a weight function between endpoints of 
strips, taking into account the direction when leaving a strip. Clearly, any feasible 
tour consists of two different matchings Mi and M2 between strip endpoints; more- 
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over, if d{Mi) and d{M2) are the total weights of the edges in the matchings, we get 
d{Mi) + d{M2) < d(T). It foUows that for an optimal matching M, we have d{M) < 
opt/2. By construction, the edges of M and the strips of 5* induce a 2-factor of the 
vertices that covers all edges. Thus, we get a cycle cover of cost at most 1.5 OPT. 

Now consider the c cycles in a cycle cover, c is at most the number of strips in 
S, which is a lower bound on the cost of an optimal tour. As the cost of merging the 
c cycles is at most 2c — 2, we get a a total cost of not more than 3.5 opt. □ 

5.2. Integral Orthogonal. As mentioned in the preliminaries, just the pixel 
count N may not be a satisfactory measure for the complexity of an algorithm, as the 
original region may be encoded more efficiently by its boundary, and a tour may be 
encoded by structuring it into a small number of pieces that have a short description. 
It is possible to use the above ideas for approximation algorithms in this extended 
framework. We describe how this can be done for the integral orthogonal case, where 
the set of pixels is bounded by n boundary edges. 

Theorem 5.6. There is a 10 -approximation of (strongly polynomial) complexity 
0{n\ogn) for computing a minimum-turn cycle cover for a region of pixels bounded 
by n integral axis-parallel segments, and a 12 -approximation of the same complexity 
for computing minimum-turn tours. In both cases, the maximum coverage of a point 
is at most 4, so the algorithms are also A- approximations on length. 

For the special case in which the boundary is connected (meaning that the region 
has no holes), the complexities drop to 0(n). 

Proof. The basic idea is to find a greedy rook cover, then use it to build an 
approximate tour. Lemma [4.121 still holds, and each strip in a star (as described in 
the previous section) will be a full strip. The approximation ratios follow as special 
cases of Theorem 15. II In this case, p = 2 and S — 4. It remains to show how we can 
find a greedy rook cover in the claimed time. 

Refer back to Figure 14.61 Subdivide the region by the n vertical chords through 
its n vertices, resulting into at most n vertical strips Xi, . . . , X„, of widths xi, . . . , Xn- 
Similarly, consider a subdivision by the n horizontal chords through the n vertices 
into at most n horizontal strips Yi , . . . , Y„ , of width yi . In total, we get a 

subdivision into at most cells Cij. Despite of this quadratic number of cells, we 
can deal with the overall problem in near-linear time: Note that both subdivisions 
can be found in time O(nlogn). For the case of a connected boundary, Chazelle's 
linear-time triangulation algorithm jl2| implies a complexity of 0(n). 

Choose any cell Cij, which is a rectangle of size Xi x yj. Then rij — min{xi,yj} 
rooks can be placed greedily along the diagonal of Cij , without causing any interfer- 
ence; such a set of rooks can be encoded as one "fat" rook, described by its leftmost 
uppermost corner (^ij , r]ij ) , and its width rij . Then the strip Xi can contain at most 
Xi — rij additional rooks, and Yj can contain at most yj — rooks. Therefore, replace 
Xi by Xi — rij , and yj by yj — rij . This changes the width of at least one of the strips 
to zero, effectively removing it from the set of strips. After at most 2n ~ 1 steps of 
this type, all horizontal or all vertical strips have been removed, implying that we 
have a maximal greedy rook cover. 

It is straightforward to see that for a fat rook at position (^^ , rjij ) and width rij , 
there is a canonical set of cycles with 10 edges each that covers every pixel that can 
be attacked from this rook. Furthermore, there is a "fat" cycle with at most 12ry — 2 
turns that is obtained by a canonical merging of the rij small cycles. Finally, it is 
straightforward to merge the fat cycles. □ 

If we are willing to invest more time for computation, we can find an optimal 
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rook cover (instead of a greedy one) . As discussed in the proof of Lemma 14.111 this 
optimal rook cover yields an optimal strip cover. An optimal strip cover can be used 
to get a 6-approximation, and the new running time is 0{n^^^ \ogn) or 0{N'^'^'^^). 

Theorem 5.7. There is an Oin?-^ log N) -time or 0{N'^-^'^^)-time algorithm that 
computes a milling tour with number of turns within 6 times the optimal, and with 
length within 4 times the optimal. 

Proof. Apply Lemma 14.111 to find an optimal strip cover of the region. (See 
Figure ^21) As described in the proof of that lemma, the cardinality of an optimal 
strip cover is equal to the cardinality of an optimal rook cover. As stated, the number 
of strips is a lower bound on the number of turns in a cycle cover or tour. 

Now any strip from u to w is covered by a "doubling" cycle with edges (u,?«), 
{w,w), (w, u), {u,u). This gives a 4-approximation to minimum-turn cycle covers. 
Finally apply CoroUarv 14 . 81 to get a 6-approximation to minimum-turn tours. 

The claim about coverage (and hence overall length) follows from the fact that 
an optimal strip cover has maximum coverage 2, hence the cycle cover has maximum 
coverage 4. □ 

By more sophisticated merging procedures, it is possible to reduce the approx- 
imation factor for tours to a figure closer to 4. Note that in the case of N being 
large compared to n, the above proof grossly overestimates the cost of merging, as 
all cycles within a fat strip allow merging at no additional cost. However, our best 
approximation algorithm achieves a factor less than 4 and uses a different strategy. 

Theorem 5.8. For an integral orthogonal polygon with n edges and N pix- 
els, there are 2. b- approximation algorithms, with running times 0{N^'^'^^ -\- n^) and 
0(n^'^ log iV -f n^), for minimum-turn cycle cover, and hence there is a polynomial- 
time Z.lb- approximation for minimum-turn tours. 

Proof. As described in Lemma 14.111 find an optimal strip cover S, in time 
q^jy2.376-| qj. O^n^-^logN). Let s be its cardinahty and let opt be the cost of an 
optimal tour, then 0PT> s. 

Now consider the end vertices of the strip cover. By construction, they are part 
of the boundary. Because any feasible tour T must encounter each pixel, and cannot 
cross the boundary, any endpoint of a strip is either crossed orthogonally, or the tour 
turns at the boundary segment. In any case, a tour must have an edge that crosses 
an end vertex orthogonally to the strip. (Note that this edge has zero length in case 
of a U-turn.) 

As in Section inl and the proof of Thcorcm l5.5l define a weight function between 
endpoints of strips, taking into account the direction when leaving a strip. Again any 
feasible tour consists of two different matchings Mi and M2 between strip endpoints, 
and for an optimal matching A/, we have d{M) < opt/2. 

Computing such a matching can be achieved as follows. Note that for N pixels, an 
optimal strip cover has 0(min{\/]V, n}) strips; by matching endpoints of neighboring 
strips within the same fat strip, we are left with 0{n) endpoints. As described in 
the proof of Lemma 14.111 the overall cost for computing the link distance between 
all pairs of endpoints can be achieved in 0(min{A^logiV, logn}). Computing a 
minimum- weight perfect matching can be achieved in time 0(max{A^^'^, ti'^}). 

The edges of M and the strips of S induce a 2-factor of the endpoints. Because 
any matching edges leaves a strip orthogonally, we get at most 2 additional turns 
at each strip for turning each 2-factor into a cycle. The total number of turns is 
2s + w{M) < 2.5-OPT. Because the strips cover the whole region, we get a feasible 
cycle cover. 
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Finally, we can use Corollary 14.81 to turn the cycle cover into a tour. By the 
corollary, this tour does not have more than 3. 75 -OPT turns. □ 

The class of examples in Example 15.91 shows that the cycle cover algorithm may 
use 2-OPT turns, and the tour algorithm may use 3-OPT turns, assuming that no special 
algorithms are used for matching and merging. Moreover, the same example shows 
that this 3.75-approximation algorithm does not give an immediate length bound on 
the resulting tour: 

Example 5.9. The class of regions shown in Figure fOI may yield a heuristic 
cycle cover with 2-OPT turns, and a heuristic tour with 3-OPT turns. 

The region consists of a "square donut" of width k. An optimal strip cover 
consists of 4fc strips; an optimal matching of strip ends yields a total of 8fc + 2 turns, 
and we get a total of 2k cycles. (In Figure IST^ a). only the vertical strips and their 
matching edges are shown to keep the drawing cleaner.) If the merging of these cycles 
is done badly (by merging cycles at crossings and not at parallel edges), it may cost 
another 4fc — 2 turns, for a total of 12fc turns. As can be seen from Figure [^Jb), 
there is a feasible tour that uses only 4fc + 2 turns. This shows that optimal tours may 
have almost all turns strictly inside of the region. Moreover, the same example shows 
that this 3.75-approximation algorithm does not give an immediate length bound on 
the resulting tour. However, we can use a local modification argument to show the 
following: 




Figure 5.2. A bad example for the 3.75-approximation algorithm: (a) Half the cycles con- 
structed by the algorithm, (b) An optimal tour. 

Theorem 5.10. For any given feasible tour (or cycle cover) of an integral or- 
thogonal region, there is a feasible tour ( or cycle cover) of equal turn number that 
covers each pixel at most four times. This implies a performance ratio of 4 on the 
total length. 

Proof. See Figure [531 Suppose there is a pixel that is covered at least five times. 
Then there is a direction (say, horizontal) in which it is covered at least three times. 
Let there be three horizontal segments (1, 1'), (2,2'), (3,3') covering the same pixel, 
as shown in Figure a) ; we denote by 1, 2, 3 the endpoints to the left of the pixel, 
and by 1', 2', 3' the endpoints to the right of the pixel. 
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Now consider the connections of these points by the rest of the tour, i.e., a second 
matching between the points 1, 2, 3, 1', 2', 3' that forms a cycle, when merged with 
the first matching (2,2'), (3,3'). This second matching is shown dashed in 

Figure lE^T b-c). We consider two cases, depending on the structure of the second 
matching. 

In the first case, there are two right endpoints that are matched, say 1' and 2'. 
Then there must be two left endpoints that are matched; because both matchings must 
form one large cycle, these cannot be 1 and 2. Without loss of generality, we may 
assume they are 1 and 3. Thus, 2 and 3' must be matched, as shown in Fieure lK^T b). 
Then we can replace (1, 1'), (2,2'), (3, 3') by (1, 2'), (2,3), (!', 3'), respectively, which 
yields a feasible tour with the same number of turns, but with some pixels being 
covered fewer times, and no pixel being covered more times than was the case in the 
original tour. 

In the other case, all right endpoints are matched with left endpoints. Clearly, 1' 
cannot be matched with 1; without loss of generality, we assume it is matched with 2, 
as shown in Figure[^3Ic). Then the cycle condition implies that the second matching 
is (1, 3'), (2, 1'), (3, 2'). This allows us to replace (1, 1'), (2, 2'), (3, 3') by (1, 3), (2, 3'), 
(1', 2'), respectively, again producing a feasible tour with the same number of turns, 
but with some pixels being covered fewer times, and no pixel being covered more 
times than was the case in the original tour. 

This can be repeated until no pixel is covered more than four times. As the above 
procedure can be carried out as part of the merging phase (i.e., after an optimal 
weighted matching has been found), the overall complexity is not affected. Further- 
more, it is straightforward to see that it also works for the case of "thick" strips, where 
N is large compared to n, by treating parallel edges in a thick strip simultaneously. □ 

5.3. Nonintegral Orthogonal Polygons. Nonintegral orthogonal polygons 
present a difficulty in that no polynomial-time algorithm is known to compute a 
minimum strip cover for such polygons. Fortunately, however, we can use the bound- 
ary tours from Section l4.ll to the approximation factor of 12 from Theorem 15.61 for 
the integral orthogonal case. 

Theorem 5.11. In nonintegral orthogonal milling of a polygonal region with 
n edges and N pixels, there is a polynomial-time A. 5 -approximation for minimum- 
turn cycle covers and 6. 25 -approximation for minimum-turn tours, with a simultane- 
ous performance guarantee of 8 on length and cover number. The running time is 
0(^2.376 ^^3)^ or 0(n2-5 log iV-^n^). 

Proof. Take the 2.5-approximate cycle cover of the integral pixels in the region 
as in Theorem 15.81 for a tour, turn it into a 3.75-approximate tour. This may leave a 
fractional portion along the boundary uncovered. See Figure l5^ 

Now add in an optimal cycle cover of the boundary which comes from Theorem 
14.41 This may only leave fractional boundary pieces uncovered that are near reflex 
vertices of the boundary, as shown in Figure Whenever this happens, there must 
be a turn of the boundary cycle cover on both side of the reflex vertex. The fractional 
patch can be covered at the cost of an extra two turns, which are charged to the two 
turns in the boundary cycles. Therefore, the modified boundary cover has a cost of at 
most 2-OPT. Compared to an optimal cycle cover of length opt, we get a cycle cover 
of length at most 4.5-OPT, as claimed. For an optimal tour of length opt, merging 
all modified boundary cycles into one cycle can be done at a cost of at most 2 turns 
per unmodified boundary cycle, i.e., for a total of i-OPT. 

Finally, the remaining two cycles can be merged at a cost of 2 turns. This yields 
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Figure 5.3. Rearranging a tour to ensure that no pixel is covered more than three times in 
each direction: (a) A set of horizontal edges that covers some pixel three times, (b) A rearranged 
tour, if the second matching between endpoints connects two left and two right endpoints. (c) A 
rearranged tour, if the second matching between endpoints connects any left with a right endpoint. 



an overall approximation factor of 3.75 + 2.5 = 6.25. The claim on the cover number 
(and thus length) follows from applying Theorem 15 . 1 01 to each of the two cycles. 

The running times follow from Theorems 14.41 and 15.81 □ 

5.4. Milling Thin Orthogonal Polygons. In this section we consider the 
special case of milling thin polygons. Again, we focus on the integral orthogonal case. 
Formally, a thin polygon is one in which no axis-aligned 2x2 square fits, implying 
that each pixel has all four of its corners on the boundary of the polygon. Intuitively, 
a polygon is thin if it is composed of a network of width- 1 corridors, where each pixel 
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fractional region along boundary 
Figure 5.4. Milling a nonintegral orthogonal polygon 

is adjacent to some part of the boundary of the region, making this related to discrete 
mining. 

5.4.1. Basics of Thin Orthogonal Polygons. Any pixel in the polygon has 
one, two, three or four neighbor pixels; we denote this number of neighbors as the 
degree of a pixel. See Figure Degree one pixels (1) are "dead ends", where 

the cutter has to make a u-turn. There are two types of degree two pixels, without 
forcing a turn (2a) or with forcing a turn {2b); in either case, applying Lemma 14.31 in 
an appropriate manner will suggest that they should be visited in a canonical way: 
after one neighbor, and before the other. Neighbors of degree three pixels (3) form 
"T" intersections that force duplication of paths. Degree four pixels (4) are the only 
pixels in thin polygons that are not boundary pixels as defined in Sectional however, 
in the absence of 2x2 squares of pixels, all their neighbors are of degree one or two. 




(1) (2a) (2b) 

— _! L 

-i r i r 

(3) (4) 
Figure 5.5. Pixel types in a thin polygon. 

In the following, we will use the ideas developed for boundary cycle covers in 
Section im to obtain cycle covers for thin polygons. The following is a straightforward 
consequence of Theorem 14.41 and Theorem 14.71 

Corollary 5.12. In thin orthogonal milling, there is an 0{n^) algorithm for 
computing a minimum-turn cycle cover, and an 0{n^) 1.5- approximation for comput- 
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ing a minimum-turn tour. 

Proof. Apply the strongly polynomial algorithm described in Theorem 14.41 for 
computing a minimum cost boundary cycle cover. By definition, this covers all pixels 
of degree one, two, and three. Moreover, degree four pixels are surrounded by pixels 
of degree one or two, implying that they are automatically covered as neighbors of 
those pixels, when applying Lemma 14.31 Using Theorem l4.7l we can turn this into a 
tour, yielding the claimed approximation factor. □ 

More interesting is that we can do much better than general merging in the case 
of thin orthogonal milling. The idea is to decompose the induced graph into a number 
of cheap cycles, and a number of paths. 

5.4.2. Milling Thin Eulerian Orthogonal Polygons. We first solve the spe- 
cial case of milling Eulerian polygons, that is, polygons that can be milled without 
retracing edges of the tour, so that each edge in the induced graph is traversed by the 
cutting tool exactly once. In an Eulerian polygon, all pixels have either two or four 
neighbors, meaning there are no odd-degree pixels. 

Although one might expect that the optimal milling is one of the possible Eulerian 
tours of the graph, in fact, this is not always true, as Example 15 .131 points out. 

Example 5.13. There exist thin grid graphs, such that no turn-minimal tour of 
the graph is an Eulerian tour. 




15 16 



5,6,7,S 5,6,7,K 
>34 ^^\9,10,11,12 /^^\ 



Figure 5.6. A thin Eulerian polygon consisting of four overlapping cycles (above). Shovin 
symbolically belovi is how to obtain an overall tour by merging the four canonical cycles: A tour 
obtained by iteratively merging cycles incurs a total of 16-1-6=22 turns (bottom left). An optimal 
tour has 16-1-4=20 turns (bottom right). 



Proof. See Figure [5l^ Observe that an optimal milling is not an Eulerian Tour. 
The best Eulerian Tour for this figure requires 22 turns, as shown symbolically in the 
bottom left of the figure: Each cycle uses 4 turns and an additional 6 turns can be used 
to connect the four cycles together. On the other hand, the optimal milling traverses 
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the edges in the internal pixel twice, both times in the same direction: The order 
of turns is 1, 2, C, 13, 16, 15, 14, D, 9, 12, 11, 10, B, 5, 8, 7, 6, A, 3, 4, 1, and the structure 
is shown symbolically in the bottom right. Thus, the optimal milling only requires 
20 turns, where each cycle uses 4 turns and an additional 4 turns connect the cycles 



By strengthening the lower bound, we can achieve the following approximation 
of an optimal tour of length opt: 

Theorem 5.14. There is an 0(n log n) (or 0{N)) algorithm that finds a tour of 
turn cost at most |-OPT. 

Proof. By applying Theorem 14. 41 we get an optimal boundary cycle cover. There 
are three observations that lead to the claimed stronger results. 

(1) For a thin polygon, extracting the coUinear strips can be performed in strongly 
polynomial time O(nlog) (or weakly polynomial time 0{N)). 

(2) For an Eulerian thin polygon, no vertices in Gp remain unmatched after 
repeatedly applying Lemma |4. 31 Instead, we get an optimal cycle cover right away. 
This cycle cover can be merged into one connecting tour by merging at pixels where 
two cycles cross each other: Let the optimal cycle cover be composed of c disjoint 
cycles, where c > 1. Let t be the cost of the optimal cycle cover. At each phase of the 
cycle-merging algorithm, two cycles are merged into one. Therefore, the algorithm 
finds a solution having cost i + 2 • (c — 1). 

(3) We can strengthen the lower bound on an optimal tour as follows. Consider 
(for c > 1) a lower bound on the cost of the optimal solution. Just like in the proof of 
Theorem 13.11 all turns in a cycle cover are forced by convex corners of the polygon, 
implying that any solution must contain these t turns. In addition, turning from one 
cycle into another incurs a crossing cost of at least one turn; thus, we get a lower 
bound of < + c. Observe that there are at least 4 turns per cycle so that t > 4c. 



5.4.3. Milling Arbitrary Thin Orthogonal Polygons. Now we consider the 
case of general thin polygons. For any odd-degree vertex, and any feasible solution, 
some edges may have to be traversed multiple times. As in Corollarv l5.12l we can 
apply Theorem 14.41 to achieve a minimum-cost cycle cover and merge them into a 
tour. Using a more refined analysis, we can use this to obtain a 4/3-approximation 
algorithm for finding an minimum-cost tour. 

Theorem 5.15. For thin orthogonal milling, we can compute a tour of turn cost 
at most I -OPT in time 0{n^ ), where OPT is the cost of an optimal tour. 

Proof. We start by describing how to merge the cycles into one connected tour. 

1. find an optimal cycle cover as provided by Theorem 14. 41 

2. repeat until there is only one cycle in the cycle cover: 

• If there are any two cycles that can be merged without any extra cost 
(by having partially overlapping coUinear edges), perform the merge. 

• Otherwise, 

— Find a vertex at which two cycles cross each other. 

— Modify the vertex to incorporate at most two additional turns, 
thereby connecting the two cycles. 

Now we analyze the performance of our algorithm. Consider the situation after 
extracting the cost zero matching edges from Lemma [4.31 This already yields a set 
K of cycles, obtained by only turning at pixels of degree two that force a turn. Let 
k denote the number of cycles in and c be the number of turns in K. Let P 
be the set of "dangling" paths at degree-one or degree-three pixels, and let p be the 



together. 



□ 



Therefore, < ^ < f- 



□ 
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number of turns in P, including the mandatory turn for each endpoint. Let M be a 
minimum matching between odd-degree vertices, and let m be the number of turns 
in M . Finally, let O be the matching between odd-degree pixels that is induced by 
an optimal tour, and let o be the number of turns in O. 

First note that P is a matching between odd-degree nodes. 

P, M and O may connect some of the cycles in K. In P a path between two 
odd-degree pixels connects the two cycles that the two nodes belong to. On the other 
hand, a path in M and O between two odd-degree nodes can encounter several cycles 
along the way, and thus it may be used to merge several cycles at no extra cost. 

Therefore, let j be the number of cycles after using P for free merging, let i be 
the number of components with P and M used for free merging, and let h be the 
number of components with P and O used for free merging. 

Note that 

(5.1) l<i<j<k, 
and 

(5.2) l<h<j<k. 

Now consider the number of cycles encountered by a path in the matching. It is 
not hard to see that this number cannot exceed the number of its turns. Therefore, 



(5.3) m>k — i>j — i 

(5.4) o> k - h> j ~ h 

If a particular matching results in x components, we would need at least another 
X turns to get a tour. Thus with O we need at least h more turns. 
Thus, for an optimal tour of cost opt, we have 

(5.5) OPT> c + p + o + h. 

Our heuristic method adds the minimum matching to C and P and merges the 
remaining components with two turns per merge, hence the cost heur of the resulting 
tour is 



(5.6) HEUR > c-l-p -I- m + 2(i - 1). 

Thus we get the following estimate for the approximation factor R> 1: 

(5.7) R < ^ + P + ^ + l^ 

c+p+o+h 

Because each cycle has at least four turns, we know that 

(5.8) c > 4fc > 4j. 

Using the fact that c + p + m + 2i > c + p + o + h (because R > 1), we see that 
the ratio on the right in (|5.7|) gets larger if we replace c in the numerator and in the 
denominator by the smaller nonnegative value 4j; thus. 
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(5.9) + + 

~ 4j+p + o + /i 

Because P is also a matching, we have 

(5.10) P>m, 
which imphcs that 

(5.11) R< 
We also know that 



Aj + 2m + 2i 
Aj + m + o + h 



(5.12) m<o. 

Using this, together with the fact that R can be assumed to be less than 2, we can 
argue that R is maximal for maximal values of m; hence, 

(5.13) R<i^':^i 

Using H5.4|) in H5.13|) . we see that R is maximal for minimal o; hence, 
f^... ^ ^ 4j + 2{j -h) + 2i 6j -2h + 2z 

Using h > o in H5.14|l and the facts that R < 2 and i < j, we get that 

6? - 2o + 2i 6j + 2i 4 
6j — o 6j 3 

□ 

The following shows that the estimate for the performance ratio is tight. 

Theorem 5.16. There is a class of examples for which the estimate of 4/3 for 
the performance ratio of the algorithm for thin orthogonal milling is tight. 

Proof. See Figure IsTI The region consists of fc = 2s + 4 cycles, all with precisely 
4 turns, s cycles without degree-three vertices, and s + 4 cycles with two degree-three 
vertices each. We get c = 4fc = 8s -I- 16 and p — 2(s + 4) = 2s + 8. Figure I^TTT bl 
shows a min-cost matching of cost m = 2(s + 4) 2s -|- 8 and one of cost o = 2s -I- 8 
that is induced by an optimum tour. As Figure \b.7i c) suggests, merging all cycles, 
odd-degree paths and matching paths is possible without requiring any further turns, 
resulting in opt = c + p + o = 12s + 32. On the other hand, using the min-cost 
matching of cost m leaves k cycles that cannot be merged for free; thus, merging two 
cycles at a time at a cost of 2 turns requires an additional cost of 2{k — 1) = 4s + 6, 
for a total of heur — c + p + m + 2{k — 1) = 16s + 38 turns, which gets arbitrarily 
close to I OPT for large s. □ 

Note that the argument of Theorem I5.1UI remains valid for this section, so the 
bounds on coverage and length approximation still apply. 
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Figure 5.7. An example with performance ratio 4/3 for our heuristic, (a) The structure of the 
example for s = 8. (b) A min-cost matching of the odd-degree vertices, and the matching induced 
by an optimal tour, (c) A portion of the optimal tour: Subtours can be merged without extra cost, 
(d) A corresponding portion of the heuristic tour: Subtours still need to be merged, which results in 
an additional cost of 2{k — 1) = 4s + 6. 



5.5. PTAS. We describe a polynomial-time approximation scheme (PTAS) for 
the problem of minimizing a weighted average of the two cost criteria: length and 
number of turns. Our technique is based on using the theory of m-guillotine subdivi- 
sions j3H|j properly extended to handle turn costs. We prove the following result: 

Theorem 5.17. Define the cost of a tour to he its length plus C times the number 
of (90-degree) turns. For any fixed e > 0, there is a (1 -\- e)- approximation algorithm, 
with running time 2'-"^'^'^ ■ N'-"^'^\ for minimizing the cost of a tour for an integral 
orthogonal polygon P with h holes and N pixels. 

Proof. Let T* be a minimum-cost tour and let m be any positive integer. Follow- 
ing the notation of pSj . we first apply the main structure theorem of that paper to 
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claim that there is an m-guiUotine subdivision, Tq, obtained from T* by adding an 
appropriate set of bridges (m-spans, which are horizontal or vertical segments) of total 
length at most ^|T*|, with length at most (1 + ^) times the length of T* . (Because 
T* may traverse a horizontal/ vertical segment twice, we consider such segments to 
have multiplicities (1 or 2), as multi-edges in a graph.) 

We note that part of Tg may lie outside the polygon P, because the m-spans 
that we add to make T* m-guillotine need not lie within P. We convert Tg into 
a new graph by subtracting those portions of each bridge that lie outside of P. In 
this way, each bridge of Tg becomes a set of segments within P; we trim each such 
segment at the first and last edge of Tg that is incident on it and call the resulting 
trimmed segments sub-bridges. (Note that a sub-bridge may be a single point if the 
corresponding segment is incident on a single edge of T* ; we can ignore such trivial 
sub-bridges.) As in the TSP method of j38j . we double the (non-trivial) sub-bridges: 
We replace each sub-bridge by a pair of coincident segments, which we "inflate" 
slightly to form a degenerate loop, such that the endpoints of the sub-bridge become 
vertices of degree four, and the endpoints of each edge incident on the interior of the 
sub-bridge become vertices of degree three (which occur in pairs). Refer to Figure l^l^ 
We let Tq denote the resulting graph. Now Tq C P, and, because T'q is obtained 
from r*, a tour, we know that the number of odd-degree vertices of Tq that lie on any 
one sub-bridge is even (the degree-three vertices along a sub-bridge come in pairs). 




I I I I I 



Figure 5.8. Definitions of sub-bridges and the graph Tq. A vertical line I is shown, which 
defines a cut. The slice of an optimal solution along I is shown, with thinner lines drawn along the 
edges of the solution (which is not intended to be an accurate instance of an optimal solution, but 
is drawn to illustrate some of the possible cases). Also shown is an enlargement of one portion of 
the cut I, showing a segment of the m-span between two portions of the boundary of P, the trimmed 
segment that forms the sub-bridge, and the portion of the resulting graph Tq in the vicinity of the 
inflated sub-bridge. 

The cost of the optimal solution T* is its length, |r*|, plus C times the number 
of its vertices. We consider the cost of Tq to be also its (Euclidean) length plus C 
times the number of its vertices. Because the number of vertices on the sub-bridges 
is proportional to their total length, and each edge multiplicity is at most two, we see 
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that the cost of T'q is 0((1 +C)/m) • |T*| greater than the optimal cost, i.e., the cost 
of T*. 

In order to avoid exponential dependence on n in our algorithm, we need to 
introduce a subdivision of P that allows us to consider the sub-bridges along an m- 
span to be grouped into a relatively small {0{h)) number of classes. We now describe 
this subdivision of P. 

By standard plane sweep with a vertical line, we partition P into rectangles, us- 
ing vertical chords, according to the vertical decomposition; see Figure ICTT al. We 
then decompose P into a set of 0{h) regions, each of which is either a "junction" or 
a "corridor." This decomposition is analogous to the corridor structure of polygons 
that has been utilized in computing shortest paths and minimum-link separators (see, 
e.g., \62\ V6'6\W^ ). with the primary difference being that we use the vertical decom- 
position into rectangles, rather than a triangulation, as the basis of the definition. 
Consider the (planar) dual graph, Q, of the vertical partition of P; the nodes of Q are 
the rectangles, and two nodes are joined by an edge if and only if they are adjacent. 
We now define a process to transform the vertical decomposition into our desired 
decomposition. First, we take any degree-1 node of Q and delete it, along with its 
incident edge; in the vertical decomposition, we remove the corresponding vertical 
chord (dual to the edge of Q that was deleted). We repeat this process, merging a 
degree-1 region with its neighbor, until there are no degree-1 nodes in Q. At this 
stage, Q has h+\ faces and all nodes are of degree 2 or more. Assume that h> 2 (the 
case ft, < 1 is easy); then, not all nodes are of degree 2, implying that there are at least 
two higher-degree nodes. Next, for each pair of adjacent degree-2 nodes, we merge 
the nodes, deleting the edge between them, and removing the corresponding vertical 
chord separating them in the decomposition. The final dual graph Q has nodes of two 
types: those that are dual to regions of degree 2, which we call the corridors, and 
those that are dual to regions of degree greater than 2, which we call the junctions. 
Each corridor is bounded by exactly two vertical chords, together with two portions of 
the boundary of P. (These two portions may, in fact, come from the same connected 
component of the boundary of P.) Each of the h bounded faces of Q contains exactly 
one of the holes of P. Refer to Figure 

Let V denote the decomposition of P just described; there is an analogous hor- 
izontal partition, H, of P into 0{h) regions. The vertical sub-bridges of a vertical 
bridge are partitioned into 0(ft) classes according to the identity of the region, r, that 
contains the sub-bridge in the vertical decomposition V. A sub-bridge intersecting a 
region t of V is called separating if it separates some pair of vertical chords on the 
boundary of r; it is trivial otherwise. (Corridor regions have only two vertical chords 
on their boundary, while junctions may have several, as many as 0(ft) in degenerate 
cases.) 

First, consider a corridor region r in V, and let a and b denote the two vertical 
chords that bound it. An important observation regarding sub-bridge classes in cor- 
ridors is this: The parity of the number of times a tour crosses a must be the same 
as the parity of the number of times a tour crosses b. The consequence is that we 
can specify the parity of the number of incidences on all separating sub-bridges of a 
given corridor class by specifying the parity of the number of incidences on a single 
sub-bridge of the class; the trivial sub-bridges always have an even parity of crossing. 

Now consider a junction region t in V. Because, in the merging process that 
defines V, we never merge a degree-2 region with a higher-degree region, we know 
that T consists of a single high-degree (> 2) rectangle, i?r, from the original vertical 
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decomposition, together with possibly many other rectangles that form a "pocket" 
attached to i?r, corresponding to a tree in the dual graph (so that removal of degree- 
1 nodes leads to a collapse of the tree and a merging of the pocket to Rr). The 
consequence of this observation is that there can be at most one separating sub-bridge 
in a junction class. (There may be several trivial sub-bridges.) 

Our algorithm applies dynamic programming to obtain a minimum-cost m-guillo- 
tine subdivision, Tq, from among all those m-guillotinc subdivisions that have the 
following additional properties: 

(1) It consists of a union of horizontal/vertical segments, having half- integral coor- 

dinates, within P. 

(2) It is connected. 

(3) It covers P, in that the center of every pixel of P is intersected by an edge of the 

subdivision. 

(4) It is a bridge-doubled m-guillotine subdivision, so that every (non-trivial) sub- 

bridge of an m-span appears twice (as a multi-edge) . 

(5) It interconnects the sub-bridges in each of a specified partition of the classes of 

sub-bridges. 

(6) It obeys a parity constraint on each of the 0{h) classes of sub-bridges: The 

number of edges of the subdivision incident on each separating sub-bridge of 
the class corresponding to a region r is even or odd, according to the specified 
parity for t. 

A subproblem in the dynamic programming algorithm is specified by a rectangle, 
having half-integral coordinates, together with boundary information associated with 
the rectangle, which specifies how the subdivision within the rectangle interacts with 
the subdivision outside the rectangle. The boundary information includes (a) 0{m) 
attachment points, where edges meet the boundary at points other than along the m- 
span; (b) the multiplicity (1 or 2) of each attachment point, and the interconnection 
pattern (if any) of adjacent attachment points along the rectangle boundary; (c) the 
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endpoints of a bridge on each side of the rectangle (from which one can deduce the 
sub-bridges); (d) interconnection requirements among the attachment points and the 
classes of sub-bridges; and (e) a parity bit for each class of sub-bridge, indicating 
whether an even or an odd number of edges should be incident to the separating 
sub-bridges of that class. There are 

^o(m) . 20ih) subproblems. At the base of the 
dynamic programming recursion are rectangles of constant size (e.g., unit squares). 

The optimization considers each possible cut (horizontal or vertical, at half- 
integral coordinates) for a given subproblem, together with all possible choices for 
the new boundary information along the cut, and minimizes the total resulting cost, 
adding the costs of the two resulting subproblems to the cost of the choices made at 
the cut (which includes the length of edges added, plus C times the number of vertices 
added) . 

Once an optimal subdivision, Tq, is computed, we can recover a valid tour from 
it by extracting an Eulerian subgraph obtained by removing a subset of the edges on 
each doubled sub-bridge. The parity conditions imply that such an Eulerian subgraph 
exists. An Eulerian tour on this subgraph is a covering tour, and its cost is at most 
0{C/m) ■ \T* \ greater than optimal. For any fixed e > 0, we set m — \C/e] , resulting 
in a (1 -l- £)-approximation algorithm with running time 0(iVO('^/'^) • 20(''). The 
techniques of I38j . which use "grid-rounded" m-guillotine subdivisions, can be applied 
to reduce the exponent on to a term, 0{C), independent of m. □ 

Remarks. It should be possible to apply a variant of our methods to obtain a 
PTAS that is polynomial in n (versus N), with a careful consideration of implicit 
encodings of tours. We note that our result relies on effectively "charging off" turn 
cost to path length, because the objective function is a linear combination of the two 
costs (turns and length). We have not yet been able to give a PTAS for minimizing 
only the number of turns in a covering tour; this remains an intriguing open problem. 

6. Conclusion. We have presented a variety of results for finding an optimal 
tour with turn cost. Many open problems remain. The most curious seems to be the 
following, which highlights the difference between turn cost and length, as well as the 
difference between a cycle cover and a 2-factor: 

Problem 6.1. What is the complexity of finding a minimum-turn cycle cover in 
a grid graph? 

This problem has been open for several years now; in fact, it is Problem #53 on 
the well-known list 17 , known as "The Open Problems Project" . Finding a minimum 
weighted turn cycle cover is known to be NP-hard for a set of points in the plane ^; 
however, the proof uses the fact that there are more than two directions for the edges 
in a convex cycle. While we tend to believe that Problem 16. II may have the answer 
"NP-complete" , a polynomial solution would immediately lead to a 1.5-approximation 
for the orthogonal case, and a (1 -I- |)-approximation for the general case. 

For various optimization problems dealing with geometric regions, there is a no- 
table difference in complexity between a region with holes, and a simple region that 
does not have any holes. (In particular, it can be decided in polynomial time whether 
a given grid graph without holes has a Hamiltonian cycle '44] , even though the com- 
plexity of the TSP on these graphs is still open.) Our NP-hardness proof makes strong 
use of holes; furthermore, the complexity of the PTAS described above is exponential 
in the number of holes. This raises the following natural question. 

Problem 6.2. Is there a polynomial-time algorithm for exactly computing a 
minimum-turn covering tour for simple orthogonal polygons? 

It may be possible to improve the performance of some of our approximation 
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algorithms. In particular, the following remains unclear. 

Problem 6.3. Is the analysis of the 3.75-approximation algorithm tight? 

We believe that it may be possible to improve the factor. We also believe that 
there is room for improvement in approximating nonintegral orthogonal milling, in 
particular by improving the cost of finding a strip cover. 

Problem 6.4. What is the complexity of computing a minimum strip cover in 
nonintegral orthogonal polygons ? 

An important tool for our approximation algorithms is a strip cover of small cost; 
finding a strip cover remains a possible approach even if strips may be parallel to 
more than two directions. This is closely related to other decomposition problems; 
see [531 fo'^ ^ survey. 

Problem 6.5. What is the complexity of computing minimum strip covers in 
nonorthogonal polygons? 

The answer may very well be "NP-hard, even for three directions" : Hassin and 
Megiddo !26! show that the problem of hitting a set of points with a minimum number 
of lines with three slopes is hard. However, their proof constructs a disconnected set 
of grid points and cannot be applied directly to milling problems. In any case, even an 
approximation would be of interest; in particular, if it achieves the following property: 

Problem 6.6. Is there a strip cover approximation algorithm for d directions 
whose performance ratio is independent of d? 

This would imply a positive result for a special case of the following, even more 
general problem. 

Problem 6.7. Can one obtain approximation algorithms for unrestricted direc- 
tions in an arbitrary polygonal domain, and an appropriately shaped cutter? 
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